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Abstract 

Short-lived  TCP  traffic  ( e.g .,  web  mice)  composes  the  majority  of  the  current  In¬ 
ternet  traffic.  Accurate  traffic  modeling  of  a  large  number  of  short-lived  TCP  flows 
is  extremely  difficult  due  to  (i)  the  interaction  between  session,  transport,  and  net¬ 
work  layers;  and  (ii)  the  explosion  of  the  size  of  state  space  when  the  number  of  flows 
is  large.  Typically,  ad-hoc  assumptions  are  required  for  the  analysis  to  be  tractable 
under  a  certain  regime. 

We  introduce  a  stochastic  model  of  a  bottleneck  ECN /RED  gateway  under  a  large 
number  of  competing  TCP  flows.  Our  main  result  shows  that  as  the  number  of  flows 
becomes  large,  the  queue  dynamics  and  the  aggregate  traffic  are  simplified  and  can 
be  accurately  described  by  simple  statistical  recursions.  These  recursions  can  be  eval¬ 
uated  independently  of  the  number  of  flows,  and  hence  the  resulting  traffic  model 
is  scalable.  Furthermore,  the  limiting  model  is  also  consistent  with  other  previously 
proposed  models  in  their  respective  regime.  Simulation  results  are  also  presented  to 
confirm  the  results. 


1  Introduction 

Due  to  the  growing  size  and  popularity  of  the  Internet,  Internet  traffic  modeling  has  become 
an  important  research  area.  Internet  traffic  consists  of  many  heterogeneous  traffic  sources, 
the  majority  of  which  utilize  TCP  congestion  control  mechanism  [4],  One  type  of  applica¬ 
tions,  such  as  FTP  and  Telnet,  are  relatively  long-lived,  while  another  type  of  applications 
are  typically  short-lived,  e.g.,  web  browsing. 

Characterization  and  modeling  of  TCP  traffic  yields  an  understanding  of  the  interac¬ 
tion  between  the  transport  layer  (TCP)  and  the  network  layer.  Such  interaction  is  well- 
understood  in  the  context  of  a  single  long-lived  TCP  flow.  However,  when  the  number  of 

‘This  material  is  based  upon  work  supported  by  the  Space  and  Naval  Warfare  Systems  Center  -  San 
Diego  under  Contract  No.  N66001-00-C-8063.  Any  opinions,  findings  and  conclusions  or  recommendations 
expressed  in  this  material  are  those  of  the  authors  and  do  not  necessarily  reflect  the  views  of  the  Space  and 
Naval  Warfare  Systems  Center  -  San  Diego. 
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TCP  flows  is  large,  straightforward  modeling  usually  results  in  a  model  that  is  not  scalable 
because  of  the  explosion  of  the  size  of  the  state  space  required  to  model  all  flows.  Fur¬ 
thermore,  for  short-lived  TCP  flows,  one  has  to  further  take  into  account  an  extra  layer  of 
dynamics,  i.e.,  the  session  layer.  These  present  two  major  obstacles  to  modeling  short-lived 
TCP  flows  accurately. 

Existing  literature  on  short-lived  TCP  traffic  modeling  usually  relies  on  ad-hoc  assump¬ 
tions,  which  render  the  model  to  be  accurate  only  in  certain  regimes.  Hollot  et  al.  model 
short-lived  TCP  flows  as  exponential  pulses,  whose  interarrival  times  are  exponentially  dis¬ 
tributed,  i.e.,  Poisson  process  [3].  They  characterize  the  statistics  of  these  exponential  pulses 
using  their  time  reversal,  which  is  described  by  the  shot  noise  process.  This  model  assumes 
that  the  short-lived  flows  last  only  a  few  round-trip  times  and  do  not  experience  packet  drops 
or  marks,  and  hence  implicitly  assumes  that  congestion  level  is  relatively  low.  Furthermore, 
flows  are  always  in  either  congestion  avoidance  (long-lived  connections)  or  slow  start  (short¬ 
lived  connections),  and  they  do  not  transition  from  one  to  the  other.  In  other  words,  they  do 
not  explicitly  model  the  session  dynamics,  where  connections  arrive  and  leave  the  network 
after  transfers  are  completed.  A  similar  approach  to  modeling  short-lived  flows  is  also  taken 
in  [7], 

On  the  other  end  of  the  spectrum,  Klierani  and  Kumar  suggest  that  as  the  bottleneck 
capacity  becomes  very  small,  the  queueing  model  for  the  bottleneck  queue  can  be  accurately 
described  as  a  processor  sharing  queue  [6].  When  the  capacity  is  large,  however,  the  processor 
sharing  model  becomes  less  accurate  because  newly  arrived  TCP  flows  cannot  fully  utilize 
their  allocated  bandwidth.  In  fact,  in  the  large  capacity  regime  these  short-lived  flows  may 
terminate  even  before  they  can  increase  their  transmission  rates  to  fully  utilize  their  allocated 
bandwidth  due  to  slow  start. 

The  shortcomings  of  these  models  suggest  a  need  for  a  unified  model  that  is  accurate  in  all 
regimes,  instead  of  being  restricted  to  one  regime  or  another.  Since  the  number  of  connections 
that  share  a  bottleneck  link  is  likely  to  be  large,  we  follow  the  approach  in  [9]  and  [10], 
which  consider  “macroscale”  modeling  of  aggregate  TCP  flows  competing  for  the  capacity  of 
a  bottleneck  link  [2],  Macroscale  TCP  models  can  be  developed  by  systematically  applying 
limit  theorems  to  derive  a  limiting  traffic  model  when  the  number  of  TCP  flows  is  large.  The 
potential  benefits  of  doing  so  are  three-fold.  First,  model  simplification  (with  the  promise  of 
scalability)  typically  occurs  when  applying  limit  theorems,  with  irrelevant  details  filtered  out 
without  relying  on  ad-hoc  assumptions.  Second,  limit  theorems  are  central  to  the  modern 
theory  of  probability,  and  as  such  have  been  the  focus  of  a  huge  literature  that  contains  a 
large  number  of  results  and  techniques.  Hence,  given  this  large  body  of  knowledge,  it  is 
reasonable  to  expect  the  existence  of  suitable  limit  theorems  (under  very  weak  assumptions) 
which  can  be  applied  to  the  situation  of  interest.  Finally,  in  the  networking  context,  resource 
allocation  problems  are  interesting  in  networks  operating  at  high  utilization,  e.g.,  when  the 
number  of  users  is  large.  In  such  a  scenario,  the  limit  behavior  will  become  increasingly  more 
accurate  as  the  number  of  users  increases. 

In  this  paper  we  extend  the  model  in  [10]  and  incorporate  an  additional  layer  of  dynamics, 
namely  the  session  layer.  We  show  that  the  queue  size  per  session  and  the  workload  per  ses¬ 
sion  brought  in  during  a  round-trip  time  converge  to  deterministic  processes  asymptotically 
as  the  number  of  sessions  increases.  Furthermore,  we  demonstrate  that  the  sessions  becomes 
asymptotically  independent,  indicating  that  the  RED  mechanism  does  alleviate  the  synchro- 
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nization  problem  among  the  connections.  The  rest  of  the  paper  is  organized  as  follows.  The 
model  and  dynamics  of  network,  transport,  and  session  layers  are  described  in  Section  2, 
followed  by  our  asymptotic  results  in  Section  3.  Section  4  gives  a  brief  discussion  on  the 
results  and  a  comparison  with  the  other  previously  proposed  models  mentioned  earlier.  A 
numerical  example  and  simulation  results  are  presented  in  Section  5  and  6,  respectively.  We 
conclude  the  paper  with  a  suggestion  for  future  work  in  Section  7. 


2  The  Model 

In  our  model,  we  have  three  layers  of  dynamics,  namely  network,  transport,  and  session 
layers,  which  interact  with  each  other  through  mechanisms  that  will  be  specified  shortly.  At 
the  lowest  level,  the  network  is  simplified  to  be  a  single  bottleneck  router  with  an  ECN/RED 
marking  mechanism  controlling  the  congestion  level.  The  traffic  injected  into  the  network 
is  controlled  by  TCP  congestion  control  mechanism  in  the  transport  layer,  which  reacts  to 
the  marks  from  the  network.  Each  TCP  connection  is  initiated  by  a  session.  A  session  can 
be  either  active  or  idle.  If  a  session  is  busy,  a  file  or  an  object  is  transferred  through  a  TCP 
connection.  A  busy  period  of  a  session  lasts  until  it  no  longer  has  any  more  data  to  transfer, 
at  which  time  it  goes  idle.  The  duration  of  an  idle  period  is  random  and  represents  the 
idle  time  between  consecutive  file  transmissions.  When  a  new  file/object  to  be  transferred 
arrives,  the  session  becomes  active  again  and  sets  up  a  new  TCP  connection.  We  now  give 
detailed  descriptions  of  the  model  for  each  layer  and  the  interaction  of  these  three  layers. 

Let  Af  =  {1,  •  •  • ,  N\  be  the  set  of  sessions  that  share  a  bottleneck  RED  gateway.  Time  is 
assumed  to  be  discrete  and  slotted  in  contiguous  time  slots  of  duration  equal  to  the  round- 
trip  delay  of  TCP  connections.  We  write  X ^  to  indicate  the  explicit  dependence  of  the 
quantity  X  on  the  number  N  of  sessions.  Equivalence  in  law  or  in  distribution  between 
random  variables  (rvs)  is  denoted  by  =st.  The  indicator  function  of  an  event  A  is  given 
by  1  [A],  and  we  use  n  (resp.  ==>„)  to  denote  convergence  in  probability  (resp.  weak 
convergence  or  convergence  in  distribution)  with  n  going  to  infinity. 

2.1  Session  Dynamics 

Each  session  i  G  Af  is  either  active  or  idle.  An  idle  session  at  the  beginning  of  time  slot 
[t,t  +  1)  does  not  have  any  packets  to  transmit  in  the  time  slot.  An  idle  session  in  time 
slot  [t,t  +  1)  becomes  active  at  the  beginning  of  time  slot  [t  +  l,t  +  2)  with  probability 
Pan  0  <  Par  <  1,  independently  of  the  past.  In  other  words,  the  duration  of  an  idle  period 
is  geometrically  distributed  with  parameter  Par  and  has  a  mean  of  1  / P„r  ■  This  attempts  to 
capture  the  dynamics  of  connection  arrivals,  where  the  interarrival  times  are  reported  to  be 
exponentially  distributed  [8].  1  Let  {£/*(£),«  G  Af]t  =  0,1,---}  be  a  collection  of  i.i.d.  rvs 
uniformly  distributed  on  [0,  1].  Let  l{p,(t+i)<-Par}  be  the  indicator  function  of  the  event  that 
a  new  file/object  arrives  in  the  time  slot  [t  +  1,£  +  2)  for  an  idle  session  i. 

Let  {Fj.(t),i  G  A f\t  =  0, 1,  -  -  -}  be  a  collection  of  i.i.d.  non-negative  integer-valued  rvs 
with  a  general  distribution  function  F.  The  workload  of  a  connection  of  session  i  that 

1  Recall  that  one  can  approximate  an  exponential  rv  A"  with  parameter  a  with  [A] ,  which  is  a  geometric 
rv  with  parameter  p  =  1  —  e~a . 
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becomes  active  at  the  beginning  of  time  slot  [t,t  +  1)  is  given  by  F,(t).  This  workload 
represents  the  total  amount  of  workload  a  TCP  connection  brings  in  before  it  is  torn  down 
rather  than  workload  brought  in  by  an  object  or  a  file.  In  other  words,  if  a  same  TCP 
connection  is  used  to  transfer  more  than  one  object  while  is  alive,  F,(t)  represents  the  total 
amount  of  workload  brought  in  by  all  objects  during  the  active  period.  We  denote  the 
remaining  workload  of  connection  i  at  the  beginning  of  time  slot  [t,t  +  1)  by  X,(t).  Clearly, 
X,(t)  =  0  if  session  i  is  idle  during  [£,£  +  1).  The  evolution  of  A", ( t )  is  given  by  the  following: 

XiN){t  +  1)  =  l{X(JV)m>0}  (X;  X  !/)  -  A\n  !/))  +  1pY!(A')p)=o}1{^h+i)<pQr}^'(^  +  1)  >  (!) 

where  A\N\t)  denotes  the  number  of  packets  injected  into  the  network  by  connection  i  at 
the  beginning  of  time  slot  [t,  t  +  1).  This  will  be  explained  in  the  following  subsection. 

2.2  TCP  Dynamics 

For  each  i  G  AC  let  w[N\t)  be  an  integer- valued  rv  that  encodes  the  congestion  window 
size  (in  packets)  at  the  beginning  of  time  slot  [t,t  +  1).  We  assume  that  the  range  of  rv 
wjN \t)  is  {0, 1,  •  •  • ,  II  max } ,  where  Wmax  is  a  finite  integer  representing  the  receiver  advertised 
window  size  of  the  TCP  connection.  We  assume  that  the  congestion  window  size  of  an  idle 
session  is  zero.  When  an  idle  session  becomes  active  at  the  beginning  of  time  slot  [t,t  +  1), 
the  congestion  window  size  of  TCP  connection  is  set  to  one  at  the  beginning  of  time  slot 
7  fi  l./  I  2).  This  models  one  round-trip  delay  for  three-way  handshake.  Here  we  describe 
how  the  congestion  window  sizes  of  active  connections  evolve.  Each  TCP  source  transmits 
as  many  of  the  remaining  data  packets  as  possible  that  are  allowed  by  its  congestion  window 
in  that  time  slot.  In  other  words,  suppose  that  connection  i  has  x[N\t)  remaining  packets 
(or  workload)  waiting  to  be  transmitted  at  the  beginning  of  time  slot  [t,t  +  l),2  the  number 
of  packets  connection  i  transmits  at  the  beginning  of  time  slot  [t,  t  +  1),  denoted  by  A^N\t), 
is  given  by 

4*7)  =  min  (wf’C),  .  (2) 

The  congestion  control  mechanism  of  TCP  operates  in  two  different  modes:  slow  start 
(SS)  and  congestion  avoidance  (CA).  A  new  TCP  connection  starts  in  SS.  In  SS,  the  con¬ 
gestion  window  size  is  doubled  every  round-trip  time  until  one  or  more  packets  are  marked. 
If  a  mark  is  received,  then  the  congestion  window  size  is  halved  and  TCP  switches  to  CA. 
The  congestion  window  size  is  limited  by  the  receiver  advertised  window  size  Umax-  Hence, 
the  evolution  of  the  congestion  window  of  connection  i  in  SS  can  be  written  as 

WiSkt  +  1)  =  min  (2U;V(/)  V  1,  Umax)  Mf](t  +  1) 

+  min  (i  _  MiiN\t+  1)):  ,  (3) 

where  a  V  b  =  max(a,  b)  and  M-N\t  +  1)  is  an  indicator  function  of  the  event  that  no  marks 
have  been  received  in  time  slot  [t,  t  +  1),  i.e. ,  M-N\t  +  1)  =  1  when  no  packet  from  Session 

2We  refer  to  a  TCP  connection  of  an  active  Session  i  by  connection  i  when  there  is  no  confusion. 
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i  is  marked  in  the  time  slot  and  M[N\t  +  1)  =  0  when  at  least  one  packet  is  marked.  The 
marking  mechanism  will  be  explained  in  more  detail  in  Subsection  2.3. 

In  CA,  the  congestion  window  size  in  the  next  time  slot  is  increased  by  1  if  no  marks  are 
received  in  time  slot  [t,t  +  1),  and  if  one  or  more  packets  are  marked  in  time  slot  [t,t  +  1) 
the  congestion  window  in  the  next  time  slot  is  reduced  by  half.  The  congestion  window  size 
in  CA  can  be  described  by  the  following: 

*0<  +  1)  =  min  +  1,  Hm„)  M[N\t  +  1) 

+  min  (1  _  Mp(t  +  l)).  (4) 

We  use  {0,  l}-valued  rvs  { .S',  x  (I ).  i  e  A/"}  to  encode  the  state  of  TCP  connections.  We 
interpret  =  0  (resp.  SlN\t)  =  1)  as  connection  i  being  in  CA  (resp.  in  SS)  at  the 

beginning  of  the  time  slot  [t,t  +  1).  Therefore,  the  complete  recursion  of  the  congestion 
window  size  can  be  written  as 

wiN>{t  +  l)  =  1(X(»)(11_H,<»)(1)>0)  (5) 

X  [S®  +  !)  +  (!-  sf  )(i))H--,c.4(i  +  1)]. 

where  the  first  indicator  function  is  used  to  reset  the  congestion  window  size  to  zero  when 
Session  i  runs  out  of  data  to  transmit  and  returns  to  its  idle  state. 

Finally,  the  evolution  of  S^N\t)  is  given  by 

SiA){t  +  1)  =  1{xlJV)  (t)~wfN)  (t)<0}  +  1{X^\t^wjN\t)>0}SiN)^MiN)^t  + 

This  equation  can  be  interpreted  as  follows.  Connection  i  is  in  SS  in  time  slot  [t  +  1 ,  t  +  2) 
if  either  (1)  there  is  no  packet  left  to  transmit  (so  the  connection  resets)  at  the  beginning  of 
the  time  slot  or  (2)  the  connection  was  active  and  in  SS  in  time  slot  [t,t  +  1)  and  received 
no  mark  in  the  time  slot.  From  (6)  we  assume  that  a  new  TCP  connection  in  SS  is  ready 
to  be  set  up  after  the  previous  connection  is  torn  down  after  finishing  its  workload,  and 
the  new  TCP  connection  becomes  active  when  a  new  file/object  arrives  initiating  three-way 
handshake. 

2.3  Network  Dynamics 

In  this  subsection  we  explain  how  packets  are  marked  to  provide  the  congestion  notification 
to  the  active  TCP  connections.  The  capacity  of  the  bottleneck  link  is  NC  packets/slot  for 
some  positive  constant  C.  The  buffer  size  is  assumed  to  be  infinite  so  that  no  packets  are 
dropped  due  to  buffer  overflow.  Thus,  congestion  control  is  achieved  solely  through  the 
random  marking  algorithm  of  the  RED  gateway. 

Let  Q  N  ( / )  denote  the  number  of  packets  queued  in  the  buffer  at  the  beginning  of  time 
slot  [t,  t  +  1).  Connection  i  injects  A^\t)  packets  into  the  network,  and  they  are  put  in  the 
buffer  at  the  beginning  of  time  slot  [t,t  +  1).  Let  the  rv 

N 

(7) 

2=1 
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denote  the  aggregate  number  of  packets  offered  to  the  network  by  the  N  sessions  at  the 
beginning  of  time  slot  [t,t+ 1).  Hence,  packets  are  available  for  transmission 

during  that  time  slot.  Since  the  bottleneck  link  has  a  capacity  of  NC  packets/time  slot, 

[Q^(f)  +  A^(£)  —  NC^+  packets  will  not  be  served  during  time  slot  [t,t  +  1),  and  will 
remain  in  the  buffer.  Hence,  their  transmission  is  deferred  to  subsequent  time  slots.  The 
number  of  packets  in  the  buffer  at  the  beginning  of  time  slot  [t  +  l,t  +.2),  Q  x  ( /  +  1),  is 
therefore  given  by 

Qm{t  +  1)  =  [Q(w)(f)  -  NC  +  A(N\t)]+  .  (8) 

Each  incoming  packet  into  the  router  in  time  slot  [£,  t  +  1)  is  marked  with  a  probability 
fC)  depending  on  the  queue  length  at  the  beginning  of  the  time  slot  [£,£+ 1).  This 

model  approximates  the  case  where  the  memory  of  the  queue  averaging  mechanism  is  long, 
which  is  the  case  for  the  recommended  parameter  settings  of  RED  [1],  We  represent  this 
possibility  by  the  {0,  l}-valued  rvs  M-j\t  +  1)  ( j  =  1, ...,  A[N\t))  with  the  interpretation 

that  M-J\t,  +  l)  =  0  (resp.  +  1)  =  1)  if  the  jth  packet  from  source  i  is  marked  (resp. 

not  marked)  in  the  RED  buffer.  To  do  so  we  introduce  the  collection  of  i.i.d.  [0, 1  ]-uniform 
rvs  {Vj.,j(t  +  1),  i,j  =  !,•••;  t  =  0, 1,  •  •  •}  that  are  assumed  to  be  independent  of  other  rvs. 
The  process  by  which  packets  are  marked  is  as  follows.  For  each  i  G  J\f  and  j  =  1,2,.. .,  we 
define  the  marking  rvs 

MS\t  +  1 )  =  1{vid(t+i  )>/(-¥)(g(JV)(b)}’ 

so  that  the  rv  M-^\t  +  1)  is  the  indicator  function  of  the  event  that  the  jth  packet  from 
source  i  is  not  marked  in  time  slot  [t,t  +  1).  The  indicator  function  of  the  event  that  no 
packets  from  connection  i  in  time  slot  [t,  t  +  1)  are  marked  can  now  be  written  as 

TWd 

Mf,(<  +  1)=  II  Wy ’(*  +  !)•  (9) 

j=  1 


3  The  Asymptotics 

The  main  result  of  the  paper  consists  of  the  asymptotics  for  the  normalized  buffer  con¬ 
tent  as  the  number  of  sessions  becomes  large.  This  result  is  discussed  under  the  following 
assumptions  (A1)-(A2): 

(Al)  There  exists  a  continuous  function  /  :  R+  — >■  [0, 1]  such  that  for  each  N  =  1,2,.. ., 

f(N\x)  =  f{N~1x)i  x>0\ 

(A2)  For  each  N  =  1,2,...,  the  dynamics  (1),  (5),  (6)  and  (8)  start  with  the  initial  conditions 
Q(N)  (0)  =  0,  wjN)  (0)  =  0,  (0)  =  1,  and  XjN)  (0)  =  0;  i  =  1, . . . ,  N.  (10) 

Assumption  (Al)  is  a  structural  condition  while  (A2)  is  made  essentially  for  technical  conve¬ 
nience  as  it  implies  that  for  each  N  and  all  t  =  0,1,---,  the  rvs  ( li'[ N * ( f ) .  S[N\t),x[N\t)), 
(u-r>«),srw,A'™  (t))  are  exchangeable.  Assumption  (A2)  can  be  omitted  but  at 
the  expense  of  a  more  cumbersome  discussion. 
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Theorem  1  Assume  that  (A1)-(A2)  hold.  Then,  for  each  N  =  1,2,...  and  t  =  0,1,..., 
there  exists  a  (non-random)  constant  q(t )  and  rvs  (W(t),X(t),  S(t ))  such  that  the  following 
holds:  (i )  The  following  convergences  take  places: 

4„,(i)  a„d  HD 

(Hf'ff),  SiA'l(l))  =!.„  X  (t),  S(i))  (12) 

(ii )  For  any  bounded  function  g  :  N3  — >■  R 

ifTS.a AKB\gmt),xm,sm)).  d3> 

Also,  jf  ElAN)(t)  4  A-  E  [min  (H-'(i),  A'(i))] .  (14) 

Moreover,  if  the  workload  distribution  F  has  a  finite  second  moment,  then 

^£,=iAf>)4,vE[A(i)]  (15) 

(in)  For  any  integer  I  =  1,2,...,  the  triplets  { (  U'/ X  (/ ).  A,  X  ( / ) .  .S';  N  (/))./  =  1 
becomes  asymptotically  independent  as  N  becomes  large,  with 

lim  P[(H',('v)(i),A,m(t),S,m(i))  =  (w„x„S,U  i=  1,  ■■-,/] 

>oo 

=  f[p[(W(t),x(t),s(t))  =  (m^iiSi)]  (16) 

2=1 

for  any  W\, . . . ,  up,  X\,  ■  ■  ■ ,  37  in  IS1  and  s  1, . . . ,  Sj  in  {0, 1}. 

In  addition,  with  initial  conditions  q( 0)  =  0,  W ( 0 )  =  0,  5(0)  =  1,  X(0)  =  0,  it  holds 


that 

q(t  +  l)  =  (q(t)-C  +  E[A(t)])+  (17) 

where  A(t )  =  min  (W(t),  X(t)).  Further,  the  recurrence 

X(t  +  1)  =st  l{g(t)=o}l-{u(t+i)<par}F(t  +  1)  +  l{x(t)>o}  (X(t)  -  A(t )) ,  (18) 

A{t)  =stmm{W{t,),X{t)) ,  (19) 

Wss(t  +  1)  =st  min  (2W(t)  V  1,  Wmax)  M(t  +  1) 

+  min^^l (I  M(t  1  1!)  ,  (20) 

WCA(t  +  1)  =st  min  (W(t)  +  1,  Wmax)  M(t  +  1) 

+  rnin^^l,lTmaxj  (1  —  M(t  +  1)),  (21) 

W(t  +  1)  =st  U  i/i  ()|  •  ( S(t)Wss{t  +  1)  +  (1  —  S(t))Wc'A{t  +  1)) ,  (22) 

S(t  +  1)  =st  l{x(t)-w(t)>a}S(t)M(t  +  1)  +  l{x(t)-m(7<o}  (23) 

holds  in  law,  where 

M(t  +  1)  =  l{V(t+i)<(i-f(q(g))A(t)j  (24) 

for  i.i.d.  [0, 1  ]-uniform  rvs  { U(t  +  1),  V(t  +  1);  t  =  0, 1, . . 

Proof:  The  proof  is  given  in  A.  ■ 
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4  Discussion 


Theorem  1  shows  that  the  dynamics  of  the  queue  at  time  t,  denoted  by  Q  ''  {!).  can  be 
approximated  by  Nq(t )  with  q(t )  determined  via  a  simple  deterministic  recursion,  which  is 
independent  of  the  number  of  sessions.  The  offered  traffic  into  the  network  during  the  time 
slot,  A(N\t) j  can  also  be  approximated  by  N  ■  E  [v4(t)] .  These  approximations  become  more 
accurate  as  the  number  of  sessions  becomes  large,  and  the  computational  complexity  does 
not  depend  on  N.  The  limiting  model  is  therefore  “scalable”  as  it  does  not  suffer  from  the 
explosion  of  state  space,  nor  does  it  require  any  ad-hoc  assumptions. 

Theorem  1  also  shows  that  the  dependency  between  each  session  becomes  negligible  under 
a  large  number  of  sessions,  i.e. ,  “RED  breaks  the  global  synchronization  when  the  number 
of  sessions  is  large.” 

Although  the  sequence  {(q(t),W(t),X(t),S(t)),  t  =  0,1,...}  is  a  time-homogeneous 
Markov  chain  with  values  in  H 1  x  }1, . . . ,  If  max}  x  N+  x  {0,1},  we  shall  not  address  here 
the  existence  of  the  steady-state  when  t  oo  as  complications  arise  due  the  fact  that  the 
first  component  is  degenerate  (i.e.,  deterministic).  However,  we  note  that  the  numerical 
calculations  for  the  limiting  model  are  very  simple.  The  number  of  steps  required  for  the 
calculation  for  each  time  step  is  independent  of  N.  We  can  determine  q(t)  through  the 
following  steps: 

(i)  Let  t  =  0.  Assume  q( 0)  =  0,  W(0)  =  0,  .Y(0)  =  0,  and  5(0)  =  1.  Use  E  [A(0)j  =  0  to 
calculate  q(l)  =  0; 

(ii)  Use  (18)  through  (23)  with  the  corresponding  q(t)  to  calculate  the  transition  prob¬ 
abilities  and  P  [(W(t  +  1)  =  X(t  +  1)  =  x,  S(t  +  1)  =  s)]  for  w  G  {0, . . . ,  H'max}, 
x  G  HNT_|_ ,  and  s  G  {0, 1}.  Then  calculate  E  [A(t  +  1)]; 

(iii)  Use  E  [A(t  +  1)]  from  (ii)  to  update  q(t  +  2)  using  (17): 

(iv)  Increase  t  by  one  then  repeat  steps  (ii)-(iv). 

We  now  consider  the  resulting  model  from  Theorem  1  in  the  regime  when  C  is  either 
very  large  or  very  small  with  the  following  assumption 

(A3)  The  marking  function  /  :  R  — >■  [0, 1]  is  monotonically  increasing  with  /( 0)  =  0  and 
lirn^oo  f(x)  =  1; 

(A4)  The  average  session  layer  input  traffic  Par  f  xclF(x)  is  strictly  non-zero. 

As  C  — >■  oo,  it  is  easy  to  see  that  q(t)  — >■  0.  Hence  the  marking  probability  per  flow 
also  converges  to  zero  from  (A3)  for  all  t.  Therefore,  each  incoming  flow  will  always  operate 
in  the  slow  start  (exponential  growth)  mode  and  hence  the  resulting  input  traffic  into  the 
network  will  become  a  multiplexing  of  (discrete-time)  Poisson  arrival  of  streams  of  random 
number  of  packets,  each  of  which  doubles  its  window  size  every  round-trip.  The  aggregate 
input  traffic  is  therefore  similar  to  the  shot-noise  processes,  which  agrees  with  [3]. 

On  the  other  hand,  as  C  — >■  0  the  queue  will  start  building  up,  and  hence  as  t.  goes  to 
oo,  q(t)  approaches  oo.  Thus,  for  large  t.  all  TCP  flows  (including  incoming  TCP  flows)  will 
experience  marking  probability  close  to  one  from  Assumption  (A3).  Therefore,  all  active 
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connections  will  be  able  to  inject  only  one  packet  per  round-trip  into  the  network  because 
every  packet  transmitted  will  be  marked,  and  TCP  congestion  window  cannot  grow  larger 
than  one.  Since  the  bottleneck  router  will  transmit  packets  non-selectively,  any  active  flow 
will  receive  roughly  equal  throughput  and  hence  the  queue  behavior  approaches  that  of 
processor-sharing,  which  is  in  agreement  with  [6]. 

5  Numerical  Example 

This  section  presents  a  numerical  example  to  study  the  behavior  of  the  queue  size  per  flow. 
The  system  and  control  parameters  are  set  as  follows: 

C  =  3  (packets/RTT),  q^in  =  2  •  N,  q*ax  =  10  •  N,pmax  =  0.1, 

The  initial  values  are  set  to  llq (0)  =  0,X2(0)  =  0,  and  5,(0)  =  1,  and  the  queue  size 
is  initialized  to  Q(0)  =  0  at  the  beginning.  The  variables  IPj(f), X^t),  and  £>*(£)  evolve 
according  to  (5),  (1),  and  (6),  respectively,  and  the  queue  size  Q(t.)  is  updated  according  to 
(8).  We  assume  the  gentle  mode  of  the  RED  gateway  in  order  to  have  a  continuous  marking 
function.  The  workload  Fi(t)  ~  Geometric(p) ,  i  =  1,  2,  •  •  •  and  t  =  0, 1,  ■  •  •,  where  p  =  0.01, 
i.e.,  E[Fi(t )]  =  100  packets.  The  idle  periods  of  sessions  are  geometrically  distributed  with 
a  mean  of  20  time  slots.  The  receiver  advertised  window  size  Wmax  is  set  to  64. 

queue  sizes  per  user 


Figure  1:  Evolution  of  queue  size  per  flow. 

Fig.  1  plots  the  evolution  of  the  queue  size  per  flow  with  N  =  100,  250,  500,  1,000,  and 
3,000  as  well  as  the  deterministic  process  q(t)  computed  through  steps  (i)-(iv)  in  Section 
4.  As  expected,  the  oscillation  in  the  queue  size  per  flow  decreases  and  the  queue  process 
converges  to  q(t )  with  increasing  N.  Although  we  have  not  formally  stated  the  convergence  of 
q(t)  to  a  steady  state,  numerical  examples  indicate  that  for  all  reasonable  set  of  parameters, 
the  queue  size  per  flow  exhibits  stationary  behavior  after  a  short  transient  period  for  all 
sufficiently  large  N. 
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6  Simulation 


In  this  section  we  verify  our  analysis  through  ns- 2  simulation  results.  In  the  simulation  we 
gradually  vary  the  number  of  sessions  from  25  to  1,000,  and  study  the  queue  behavior.  Other 
system  parameters  are  scaled  with  the  number  of  sessions  N  as  follows: 

CN  =  0.24  •  N  Mbps,  BN  =  25  •  N, 

Qmin  =  2  •  iV,  qfnax  =  10  -N,  pmax  =  0.1  , 

where  CN  is  the  bottleneck  link  capacity,  qffin  and  qfnax  are  the  threshold  values  used  by  the 
RED  gateway  [2],  and  is  the  buffer  size.  The  receiver  advertised  window  size  is  set  to  64 
packets,  and  the  packet  size  is  set  to  1,000  bytes.  The  exponential  averaging  weight  of  the 
RED  gateway  is  set  to  0.02/iV  in  order  to  have  a  similar  time  constant  in  all  cases.  A  session 
generates  a  workload  that  is  exponentially  distributed  with  a  mean  of  100  packets,  and  the 
interarrival  times  of  the  new  workloads  for  each  connection  are  exponentially  distributed 
with  a  mean  of  3.3  seconds.  When  a  session  runs  out  of  data  to  transfer,  it  terminates  the 
TCP  connection.  A  new  TCP  connection  is  initiated  by  the  session  when  the  next  workload 
arrives  for  the  session.  The  round-trip  propagation  delays  of  the  sessions  are  randomly 
selected  from  [52,  121.5]  ms,  with  a  mean  of  approximately  87  ms.  The  gateway  implements 
the  RED  mechanism  with  ECN  option,  and  the  gentle  mode  of  the  RED  mechanism  is  turned 
on  so  that  the  marking  function  is  continuous  in  the  average  queue  size.  Also,  the  drop -front _ 
option  of  the  RED  gateway  is  enabled,  i.e.,  the  RED  gateway  marks  the  packet  at  the  front 
of  the  queue  rather  than  the  packet  that  has  just  arrived,  in  order  to  reduce  the  feedback 
delay.  The  TCP  connections  are  Reno  connections  that  are  ECN  capable,  i.e.,  connections 
react  to  ECN  marks  set  by  the  gateway. 


queue  size  per  user 


time  (sec) 


Figure  2:  Evolution  of  queue  size. 

Fig.  2  shows  the  evolution  of  the  queue  size  per  flow  with  different  number  of  sessions.  As 
one  can  see  the  magnitude  of  oscillation  in  the  queue  size  per  flow  decreases  with  increasing 
N,  and  the  queue  size  per  flow  converges  to  a  deterministic  process  as  predicted  by  our 
results. 
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7  Conclusions 


In  this  paper,  we  have  developed  a  stochastic  model  for  general  TCP  flows,  which  has  taken 
into  account  the  interaction  between  three  layers,  namely  the  network,  transport,  and  session 
layers.  The  resulting  model  is  scalable  and  is  more  accurate  as  the  number  of  sessions  grows 
large.  A  sharper  approximation  can  also  be  developed  with  a  central  limit  theorem-type 
complement  similar  to  [10]. 

While  the  limiting  result  applies  only  to  TCP  flows  with  identical  round-trip,  it  will  be 
of  use  in  a  number  situations,  e.g.,  the  buffer  dimensioning  problem  in  an  intercontinental 
Internet  link,  where  the  intercontinental  link  is  typically  a  bottleneck,  its  large  propagation 
delay  dominates  the  round-trip  delays  of  the  connections,  and  the  number  of  sessions  is 
extremely  large.  Also,  our  limited  simulation  results  indicate  a  similar  limiting  behavior  for 
connections  with  heterogeneous  round-trip  delays. 

Although  we  have  yet  to  prove  the  existence  of  a  steady  state  regime  for  the  limiting  re¬ 
cursion  identified  here,  the  limited  simulation  results  suggest  the  existence  of  such  a  steady 
state  under  some  conditions.  Future  work  on  this  class  of  models  also  includes  the  incorpo¬ 
ration  of  random  round-trip  delays  and  non-responsive  UDP  flows. 
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A  A.  Proof  of  Theorem  1 

A.l  Some  simple  and  useful  facts 

To  facilitate  the  presentation  of  the  proof  of  Theorem  1,  we  begin  with  a  few  simple  and 
useful  facts.  Fix  i  =  1 , . . . ,  Ar  and  consider  an  arbitrary  bounded  mapping  g  :  IS3  — >■  R: 
Through  careful  case  analysis,  it  follows  from  (1),  (5)  and  (6)  that 

9 ( (f  +  1 ) ,  X^  (t  +  1 ) ,  SjN^  (f  +  1 ) ) 

=  1{X[JV)(P  =  0}^  (0’  1{Ui(t+l)<Par}Fi(t  +  !))  t) 

+  1{x]JV)(t)>uAv)(p,AA')(p=o} 

x  (MjN\t  +  l)Fi  (w!N)(t),xlN)(tj)  +  (l  -  MjN)(t  +  lj)  F2  (\\fT\t),x\N\t))) 

+  h  vT  v)(0Ap  :t)  1 1 

x  +  1)F3  +  (l  -  MiiN)(t+  1))  F2  (w!*\t),xlN\tj)) 

+  1  { n<  .V]  v  :/:■  U;  v  (^5) 

where  the  IN  x  IN  mappings  iy] .  iyj  and  Fg  are  associated  with  g  and  defined  as  follows: 

F;  =  g  (min  +  1.  .  Af°(«)  -  tlf>(t),o) 

f  >(i))  =  g  (min  ] ,  M,’max  j 

F?  (wfAlXR >(*))  =  g  (min  (2W,:(S ' \t)  V  1,  WmK)  ,  X<N,(t)  -  W<N,(t ),  l)  .  (26) 

Despite  the  complicate  form  of  (25),  it  can  be  constructed  intuitively.  First,  note  that  the 
state  of  the  triplet  (xjN\t  +  1),  W-N\t+  1),  5?-/V)  (t  +  1))  depends  on  the  state  of  the  triplet 
at  time  t.  For  example,  if  xj N \t)  =  0,  then  W-N\t  +  1)  has  to  be  zero  and  S'f N 1  ( t  +  1)  is 
always  one  while  the  value  of  x[N\t  +  1)  depends  on  whether  a  new  request  arrives  or  not, 
hence  we  have  the  first  term  on  the  RHS  of  (25).  The  rest  of  (25)  can  also  be  understood 
in  a  similar  fashion  and  the  detail  will  be  omitted  here.  This  simple  fact  that  any  arbitrary 
mapping  of  the  triplet  at  time  t  +  1  can  be  expanded  as  a  function  of  the  triplet  at  time 
t  with  additional  random  elements  that  are  independent  of  any  events  up  to  and  including 
time  t  is  crucial  later  in  the  proof  of  Theorem  1. 

Let  Tt  denote  the  cr-field  generated  by  the  rvs  {Q^(0),  1L;  x  (0).  xf^O),  £y(s),  F*(s), 
Vi{s),Vid{s)  ,  i,j  =  1,2,...  :  s  =  }.  with  the  rvs  Q(N]{t),  xjN)(t),  SjA)(t)  and 

W-N\t)  (i  =  1  being  all  ^-measurable,  it  holds  under  the  enforced  independence 

assumptions  that 

E  I  1)A  1  fiN)(QiN)(t)),  j  1.2.... 
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so  that 

E  [iWf>(t  +  l)|7-t]  =  z\N\t) 

(27) 

by  conditional  independence,  where  we  have  set 

AN) ,  v 

z\N\t)  =  {l  -fiN\QiN)(t)))  < 

(28) 

It  is  now  clear  that 

Mi  \f  +  !)  =*t  1{v;:p+i)<4JV)p)}' 

(29) 

It  readily  follows  from  (25)  that 


E  [9(H-f '(i  +  1),  Xf\t  +  1),  Sfi(t  +  1))|^,] 

=  1{xlN)(t)= 0}  eK° 

+  Vr  >«>»-<»><  ,i,sr>«)=o)  +  (i 

+  1{x<»>(o>h-<"'(,i4”'«).1)  {z(tN)(t)F^w»\t),xr\t))  +  (l 

=  F,(Z\K\t),  Wf  >(f),  -Yf>(i),  Sf’fi)) 


(30) 


where  the  mapping  Fg  :  [0, 1]  x  IN  x  IN  x  {0.1}  — >■  R  is  associated  with  g  through 

Fg{z,  w,  x,  s)  =  l{z=o}E[#(o,  l{cri(i+i)<Por}i?i(i  +  l),l)] 

+  Mx>w}Ms=o}  {zFg{w,x)  +  (1  -  z)Fj}(w,x)'j 
+  l{.c>w}l{s=i}  ( zF*{w ,  x)  +  (1  -  z)Fg(w ,  re)) 

+  l{o<.c<w}5,(0,  0, 1).  (31) 

We  note  that  E  \cj  (0,  l{[/,:p+i )<par}Fi{t  +  1),  l)]  always  exists  and  is  finite  because  the  map¬ 
ping  g  is  bounded.  Further,  the  mapping  Fg  is  continuous  with  respect  to  the  product 
topology  on  [0, 1]  x  IN  x  IN  x  {0, 1}. 

Upon  taking  expectations  on  both  sides  of  (30)  we  see  that 


E  [g  (H'fV  +  l)Wf ’(«  +  1).  sfU  +  1))]  =  E  [i?  (Zg’li),  wlN\t),  X\K '(«),  5f : '(«))]  -(32) 


A. 2  A  Weak  Law  of  Large  Numbers 

We  introduce  the  following  terminology  to  facilitate  the  discussion:  For  each  t  =  0,1,...,  the 
statements  [A:t],  [B:t],  [C:t]  and  [D:t]  below  refer  to  the  following  convergence  statements: 

[A:t]  For  some  non-random  q(t ),  it  holds  that 

~ t  ;.\(/(/):  (33) 
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[B:t]  For  some  {1, . . . ,  W,rmax) -valued  rv  W(t),  non-negative  integer-valued  rv  X(t),  and 
{0,  l}-valued  rv  S(t),  it  holds  that 

(H'-1<A',(<),Atv,(i).S‘A')(())  =s-jv  (W(t),X(t),  S(t)) ;  (34) 

[C:t]  For  any  integer  1  =  1.2,...,  the  rvs  j  ( 1 1  (  t  j .  \"  ’  i  t ),  .S’)  V 1 !  /  j  j  ,  /  I  , . . . .  /  }  become 

asymptotically  independent  with  large  N  as  described  by  (16)  where  (W(t),X(t),  S(t )) 
are  the  rvs  occurring  in  [B:t]; 

[D:t]  For  any  bounded  mapping  g  :  IN'3  — >■  R,  the  convergence  (13),  (14)  hold  with  {W(t), 
X (£),  S(t))  the  rvs  occurring  in  [B:t].  Moreover,  if  the  file  arrival  distribution  has  a 
finite  second  moment,  then  the  convergence  (15)  also  holds. 

With  the  help  of  a  series  of  lemmas,  we  shall  prove  the  validity  of  the  statements  [A:t]- 
[D:t]  for  all  t  =  0, 1, . . ..  We  do  so  by  induction  on  t  and  in  the  process  we  establish  Theorem 

1. 

Lemma  1  Under  (Al),  if  [A:t]  and  [B:t]  hold  for  some  t  =  0, 1, . . then  [B:t+1]  holds 
with  W(t  +  1),  X(t  +  1),<5(£  +  1)  related  in  distribution  to  W(t),  X(t),  S(t)  by  (22),  (18), 
and  (23). 

Proof:  Together  the  convergence  [A:t]  and  [B:t]  imply  [5,  Thm.  5.28,  p.  150]  the  joint 

convergence  (N~1Q(-N')(t),wlN\t),x[N\t),S{N\t))  =>jv  W(t),  X(t),  S(t)).  Next  the 
continuity  of  the  mapping  /  implies  that  of  ( y,  w,x )  — >■  (1  —  f(y))mm(w'x)  0n  R+  x  [0,  oo)  x 
[0,  oo),  so  that 


(zSK)(<),^1W(<),.Yf)(<),sS®(<))  (Z(t),W(t),X(t),S(t))  (35) 

by  the  Continuous  Mapping  Theorem  [5,  Thm.  5.29,  p.  150]  with 

Z(t)  =  (1  -  f(q(t))r^w^xm  . 

Consider  (32)  for  any  bounded  arbitrary  mapping  g  :  IS3  — >■  R,  and  recall  that  the 
mapping  Fg  defined  by  (31)  is  continuous  on  [0,1]  x  Is1  x  IN  x  {0,1}.  Consequently,  the 
Continuous  Mapping  Theorem  can  again  be  invoked  to  yield 

F’a(4A',(<)-lh'V,C).A'!'v,(<).s!-v,(i))  Fa(Z(<),H-(i),A(i),S(i)),  (36) 

whence 

lim  EklzWliMFf'liJ.A'riO.Sr'li))]  =E[Fa(Z(i),FF(i),A'(i),S(«))]  (37) 

JX —y  OO  L  J 

by  the  Bounded  Convergence  Theorem  [5,  Thm.  4.16,  p.  108].  Combining  (32)  and  (37)  we 
get 

lim  E  [g(\\f](l  +  1  ),X?)(t  +  1),  S[N\t))}  =  E  [Fg(Z(t),  1F(/)..V(/).  S(t))[  (38) 
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and  since  the  bounded  mapping  g  is  arbitrary,  it  follows  immediately  that 

( w[N\t  +  1),  x[N\t  +  1),  s[N\t  +  1))  =>$  (W(t  +  1),  X(t  +  1),  S(t  +  1)) 

for  some  {1, . . . ,  il'Tnax}  valued  rv  W(t  +  1),  non-negative  integer-valued  X(t  +  1)  and  {0,1}- 
valued  rv  S(t  +  1)  with 

E  [g  ( W(t  +  1),X (t  +  1),  S(t  +  1))]  =  E  [Fg  ( Z(t ),  W(t),X(t),  S(t))]  .  (39) 

A  moment  of  reflection  and  a  comparison  to  the  analysis  in  (30)-(32)  will  convince  the  reader 
that  (39)  is  equivalent  to  (18)-(23).  ■ 


Lemma  2  Under  ( A1 ),  if  [A:t]  and  [D:t]  hold  for  some  t  =  0,1,...,  then  [A:t+1]  also 
holds. 


Proof:  From  [A:t]  and  [D:t]  (  specifically,  (14)  ),  we  conclude  that 

''  A  Si  i  •  '!  \t)  “A  n  q{t)  ~  C  +  E  [A(t)]  (40) 

and  the  desired  result  is  a  simple  consequence  of  the  continuity  of  the  function  x  — >■  x+  as 
we  note  that  since 


QiN)(t  +  l) 
N 


Q(N)(t) 


N 


-C  +  ^ElA\t) 


for  all  N  =  1,2, _ 

The  proof  of  Lemma  2  also  shows  that 

—  ~ Nq(t  +  1) 

with  non-random  q(t  +  1)  determined  by  (17). 


Lemma  3  Under  ( A1)-(A2 ),  if  [A:t],  [B:t]  and  [C:t]  hold  for  some  t  =  0,1,...,  then 
[C:t+1]  also  holds. 

Proof:  We  first  observe  that  for  a  fixed  N,  the  triplets  (WjN\t),  x!-N\t),  Sjw\t)  are 

coupled  only  through  the  marking  probability  which  depends  only  on  Q  ''  {!).  Fix  a  positive 
integer  I.  The  rvs  V\ (t  +  1), . . . ,  \ /(t  +  1)  are  i.i.d.  [0,  l]-uniform  rvs  which  are  independent 
of  T\.  Thus,  upon  making  use  of  the  representation  (3)-(4)  with  (29),  we  see  that  the  rvs 
( w[N \t.  + 1 ) ,  x[N \t,  + 1 ) ,  s[N \t  + 1 )),...,  ( w\N \t.  + 1 ) ,  x\N \t  + 1 ) ,  S[N \t.  + 1 ) )  are  mutually 
independent  given  T\.  Consequently,  for  arbitrary  bounded  mappings  gi, . . .  ,gj  :  N  — )■  R, 
we  get 


E 


n  9.(»a)(* + a,  x?\t + a,  5f  y +i)i^, 


4  N) 


(N), 


.2=1 


=  II E  [9.(W',<")(*  +  1).  A'f  >((  +  1),  sf  >(«  +  1)1^,] 

2=1 

=  n  Fv.UUm.  + a-  sf’t* +i)) 

2=1 
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with  the  help  of  (30)  and  (31). 

Now  it  follows  from  (16)  in  [C:t]  that  the  joint  convergence 

(Kf  >((),  A-rV),  sf '(0  ....  Hf '((),  Af ’>(*),  S‘K'(*)) 

=S-jv  (W'i(<),  A St A',(<),  S,(<)) 

holds  with  limiting  rvs  (Wi(Z),  A^(t),  Si(t)) . . . ,  (W/(t),  Xj(t),  5/(2 ))  which  are  i.i.d.  rvs  each 
distributed  according  to  (W (t) ,  X (t) ,  S (t)) .  As  in  the  proof  of  Lemma  1,  the  arguments 
leading  to  the  convergence  (36)  also  lead  to 

(Fsl(z[N)m,  >(*),  a'!a'>(().  s<a'V)).  ....  M-;  ">(*),  sf>m 

=>«(Fg,[ZM,Wl{t),Xl{t),SI[t)) . F„(Z,  (t),W,(t)),Xr(t),Sr(t)) 

where  the  limitingrvs  W\(t),  _Yi(2),  S'i(t)), . . . ,  (Zi(t),  Wi(t),Xi(t),  5/(2 ))  are  i.i.d.  rvs 

each  distributed  according  to  the  pair  (Z(t),W(t),X(t),S(t)).  Therefore,  by  the  Bounded 
Convergence  Theorem, 


lim  E 

W->oo 


=  lim  E 

N— >-oo 


[|  9.(»f ’(«  +  1),  ATf  +  1),  Sf  >(f  +  1)) 

i— 1 

II  FS,  (#’(«)  ■ 1  l'f °  ( t ) ,  A,1' (t ) ,  Sf 1  ( f ) ) 

2=1 

2=1 

2=1 

=  n  E  \9i{Wi{t  +  1),  Xi(t  +  1),  5,(2  +  1))] 


=  E 


(41) 


2=1 


where  the  last  equality  made  use  of  the  relation  (39).  The  desired  result  [C:t+1]  now  follows 
from  (41)  given  that  the  mappings  gi, ...  ,gT  are  arbitrary.  ■ 

Lemma  4  Under  (A1  )-(A2),  if  [A:t] ,  [B:t]  and  [C:t]  hold  for  some  t  =  0, 1, . . .,  then  [D:t] 
holds. 

Proof:  Pick  a  mapping  g  :  N3  — >■  R.  We  begin  by  observing  that  under  (A2)  the  rvs 

(WfN\t),xjN\t),  SjN\t))*f  i  =  1, . . . ,  N  are  exchangeable.  As  a  result,  we  get 


™r[^E£i»Wm(i).Al‘v)(<),s;'')(i)) 

N 

=  N -2  X  var[9(Hf '>(<),  Af°(«),  Sf ’(«))] 


2=1 


N 


+  N~2  E  cov[g(Wr\t),XlN\t),SlN\t)),g(W^\t),X]N\t),Sf\t))] 

=  JV-1Tar[<,(Hfv)(<),A'|A',((),S|'v)(<))] 

+  ^-cov[9(Hf '(«),  Af°(i),  s!N)(i)).«(  Wl"1 ’  W,  aT°(<).  Sf ’(*))].  (42) 
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Now  let  N  go  to  infinity  in  (42):  The  validity  of  [C:t]  and  the  Bounded  Convergence 
Theorem  already  imply 


Jim  cov[9(H'f  >(<),  X ST'(i))] 

N— s-oo 

=  co  Y[g(W1(t),X1(t),S1(t)),g(W2(t),X2(t),S2(t))]  =  0 
by  asymptotic  independence.  On  the  other  hand, 

limsupvar[g(Wf^(£),  A^^),  Sf^(£))]  <  oo 

N— s-oo 

since  g  is  bounded. 

Combining  these  observations  we  readily  see  that 


AN)r+\  \  (v !  /  /  \  c(N), 


lim  var  £*=i g{W}N){t),xlN){t),SlN\t)j\  =  0, 

AT— ino  L  J 


whence,  by  Chebysliev’s  Inequality, 


4E.=rf(CN!(<)Wf)(i),sf)(i))  -  E  [iE£1«(»Tim(i),.Yf®(<),sf)(t))]  4  *0. 


This  last  convergence  is  equivalent  to 


bv  exchangeability,  and  the  desired  convergence  result  (13)  is  now  immediate  once  we  remark 
under  [B:t]  that  linqv_!.00  E  [g{WiN\t),  x[N\t),  S[N\t))  =  E  [g(W(t),  X(t),  S(t))].  It  is 

then  straightforward  to  get  (14)  with  g(W^N\t),  X^N\t))  =  A^\t)  =  mm(W,jA\t),  X^\t)) 
which  is  bounded  by  IT  max- 

Finally,  if  the  file  arrival  distribution  F  has  a  finite  second  moment,  the  convergence  (15) 
follows  from  the  dominated  convergence  theorem  on  (43)  with  g(W-A \t),  xjA\t,):  slN\t))  = 
xjA\t)  and  note  that  xjA\t)  is  dominated  by  a  random  variable  with  the  distribution  F 
which  has  a  finite  second  moment  .  ■ 

We  now  conclude  with  a  proof  of  Theorem  1:  We  first  note  that  under  (A1)-(A2)  the 
statements  [A:t]-[D:t]  trivially  hold  for  t  =  0.  Moreover,  if  [A:t]-[C:t]  hold  for  some 
/  0,  1 . . . .,  then  so  do  the  statements  [D:t]  [B:t+1],  [A:t+1]  and  [C:t+1]  by  Lemma  4, 

Lemma  1,  Lemma  2  and  Lemma  3,  respectively.  Consequently,  the  statements  [A:t]-[D:t] 
do  hold  for  all  t  =  0, 1, . . .  by  induction  and  the  validity  of  Claims  (i)-(iii)  of  Theorem  1  is 
now  established.  The  dynamics  (17)  is  a  byproduct  of  the  proof  of  Lemma  2,  while  (18)-(24) 
are  already  contained  in  Lemma  1. 
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